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Abstract. We review our work on the application of the renormalization-group method to obtain first- 
and second-order relativistic hydrodynamics of the relativistic Boltzmann equation (RBE) as a dynamical 
system, with some corrections and new unpublished results. For the first-order equation, we explicitly obtain 
the distribution function in the asymptotic regime as the invariant manifold of the dynamical system, which 
turns out to be nothing but the matching condition defining the energy frame, i.e., the Landau-Lifshitz 
one. It is argued that the frame on which the flow of the relativistic hydrodynamic equation is defined 
must be the energy frame, if the dynamics should be consistent with the underlying RBE. A sketch is also 
given for derivation of the second-order hydrodynamic equation, i.e., extended thermodynamics, which is 
accomplished by extending the invariant manifold so that it is spanned by excited modes as well as the 
zero modes (hydrodynamic modes) of the linearized collision operator. On the basis of thus constructed 
resummed distribution function, we propose a novel ansatz for the functional form to be used in Grad 
moment method; it is shown that our theory gives the same expressions for the transport coefficients as 
those given in the Chapman-Enskog theory as well as the novel expressions for the relaxation times and 
lengths allowing natural interpretation. 
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1 Introduction 

The dynamical evolution of the hot and/or dense QCD 
matter produced in the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory can be well 
described by relativistic hydrodynamic simulations [1,2]. 
It seems to be the case also for the created matter in 
LHC (Large Hadron Collider) in European Organization 
for Nuclear Research (CERN); see, for example, [3,4]. The 
suggestion that the created matter at RHIC may have 
only a tiny viscosity prompted an interest in the origin of 
the viscosity in the created matter to be described using 
the relativistic quantum field theory and also the dissi- 
pative hydrodynamic equations. We note that since the 
created matter expands, the proper dynamics for the de- 
scription may change from hydrodynamics to kinetic one 
and vice versa [4-9] . The hydrodynamics is also relevant to 
the soft-mode dynamics [10-12] around the possible crit- 
ical point(s) in QCD phase diagram [13]; see [14] for the 
latest up date. 

However, the theory of relativistic hydrodynamics for 
viscous fluids is still under debate. In fact, we can indicate 
the following problems: (1) There are ambiguities in the 
definition of the flow velocity [15-18]. (2) In the Eckart 
(particle) frame, there arises an unphysical instabilities of 
the equilibrium state [19]. (3) The so called first-order 



equations lack in causality, i.e., some components of the 
hydrodynamic equations are of parabolic nature [20-23] . 

Taking the relativistic Boltzmann equation (RBE) [20, 
21] as a typical kinetic equation, we have been exploring 
the basic problems with the relativistic hydrodynamics 
[24-26] . We note that such an approach is important also 
for a systematic analysis of RHIC /LHC data, because the 
proper dynamics for the description may change from hy- 
drodynamics to kinetic one and vice versa, as mentioned 
above. 

It is conjectured [27, 28] that the non-equilibrium pro- 
cess evolves through some stages of hierarchical dynamics: 
In the beginning of the time evolution of an isolated pre- 
pared state, the whole dynamical evolution of the system 
will be governed by Hamiltonian dynamics that is time- 
reversal invariant. As the system gets old, the dynamics is 
relaxed into the kinetic regime where the time-evolution 
system is well described by kinetic equations which de- 
scribe a coarse-grained slower dynamics: The Boltzmann 
equation for the one-body distribution function is one of 
them [28] . Usually the original time- reversal invariance is 
lost in the description by the kinetic equation through the 
coarse-graining. As the system is further relaxed, the time 
evolution will be described in terms of the hydrodynamic 
quantities, i.e., the flow velocity, particle-number density, 
and local temperature. In this sense, the hydrodynamics is 
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the far-infrared asymptotic dynamics of the kinetic equa- 
tion. 

Thus, for obtaining the proper relativistic hydrody- 
namic equation, it is a legitimate and natural way to start 
with the RBE which is Lorentz invariant and expected 
to be free from causality problem [20,21]; moreover, ap- 
parent instability is not known for numerical simulations 
of the RBE, as far as we are aware of, and the stability 
is proved at least for the linearized version of it [29,30]. 
For analyzing the problems (1) and (2) first, we derive 
hydrodynamic equation [24,26] from the RBE. We note 
that the problem is a typical reduction problem of a dy- 
namical system in the far-infrared long- wave length limit. 
So we need a powerful reduction theory for our purpose, 
and we shall take the renormalization-group (RG) method 
[31-33] as such a powerful one. The reduction of dynamics 
can be viewed as a construction of an invariant /attractive 
manifold [34,35], and it has been shown [32,33] that the 
RG method can be nicely formulated as an elementary 
method for constructing the invariant manifold of a given 
dynamical system. 

In this article, we also report on our attempt [36] to 
examine the causality problem (3) by deriving the so called 
extended thermodynamics [38-40] : Namely, we derive meso- 
scopic dynamics of the RBE by constructing the invari- 
ant/attractive manifold incorporating some fast modes as 
well as the zero modes of the linearized collision operator. 
It turns out that our theory leads to the same expressions 
for the transport coefficients as given by the Chapman- 
Enskog method [41] and also novel formulas of the relax- 
ation times in terms of relaxation functions, which allow 
a natural physical interpretation of the relaxation times. 
Moreover, the distribution function which is explicitly con- 
structed in our theory provides a new ansatz for the func- 
tional form of the distribution function in the Grad theory 
[42]. 

2 Introduction to renormalization-group 
method by an example 

Our approach is heavily based on the reduction theory of 
dynamics called the renormalization-group (RG) method 
[31-33], and the reliability of our theory is assured by that 
of the method. It is nice [33] that the RG method can be 
formulated as an elementary way of construction of the 
invariant/attractive manifold of dynamical systems; it not 
only leads to asymptotic dynamics of a given equation but 
also extract explicitly the differential equations governing 
the would-be constants appearing in the solution to the 
differential equation. In this section, we make an account 
of the RG method using a simple non-linear equation. 
Let us take the Van der Pol equation as an example: 

x + x = e(l -x 2 )x, (2.1) 

with e being small. 

Let x(t; t ) be a local solution around t ~ Vio, and 
represent it as a perturbation series; 

x(t; t ) = x (t; t Q ) + ex 1 (t; t ) + e 2 x 2 (t; to) + ■■■{2.2) 



In the RG method, the initial value W(t ) is to constitute 
the desired (approximate) solution in a global domain and 
make the invariant manifold of the system. We suppose 
that an exact solution is given by x(t) and the initial value 
of x(t; to) at t = to is set up to be x(to), i.e., 

W(t ) = x(t ; to) = x(t ). (2.3) 

The initial value as the exact solution should also be ex- 
panded as 

W(to) = Wo(to) + eW 1 (to)+e 2 W 2 (to) + ---. (2.4) 
The zeroth-order equation reads 

Cxo = 0, (2.5) 

,2 

with Co = j[£z + 1- The solution may be expressed as 

x (t; to) = A(t ) cos(t + 9(to)), (2.6) 

with the initial value W(to) — xo(to; to) = A(to) cos(io + 
6 (to)). Note that the integral constants A and 9 may de- 
pend on the initial time to- The integration constants 
A(t ) and 9(to) will parametrize the global solution and 
correspond to the hydrodynamic variables which parametrize 
the distribution function in the local equilibrium. 
The equation for xi reads 

/ A 2 \ A 3 
Cx 1 = -A[\- — \ sin 0(f) + — sin 3</>(f), (2.7) 

with (j)(t) = t + 6 (to)- Notice that the first term in r.h.s. 
is a zero mode of C, and hence the special solution to this 
equation contains a secular term that is expressed as a 
product of t and a zero mode of C Since we have sup- 
posed that the initial value at t = to is on an exact solu- 
tion, we should make the corrections from the zeroth-order 
solution as small as possible. This condition is realized by 
letting the secular terms vanish at t — to, which is possible 
because we can freely add zero mode solutions to a special 
solution. Thus, we have the first-order solution as 

A ( A 2 \ 
ii(t;t ) = (t-to)- ( 1-— J sin 0(f) 

A 3 

~ 32 Sin30(<) ' (2 ' 8) 

with the initial value at t = t , Wi(t ) = xi(t ;t ) = 
-A 3 (to)/'A2 ■ sin30(f o ). 

If we stop here, we have the perturbative solution; x, = 
xo + exi, which should be valid in a local domain t ~ to 
but becomes invalid in the global domain where \t — t \ 
can be large, due to the secular term. 

We shall now take a geometrical point of view [32]: 
The function x(t;t ) corresponds to a curve drawn in the 
(t, x) plane for each to; in other words, we have a family of 
curves represented by x(t; to) in the (t, x) plane; a member 
of the family is parametrized by to, and each member is 
close to an exact solution in the neighborhood of t = to- 
Thus, an idea is that the envelope curve of the family of 
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curves should give a global solution. The classical theory of 
envelope curve says that the envelope can be constructed 
by solving the following equation, 



dx 
dt 



0. 



t =t 



(2.9) 



which is called the RG equation [31]; we have here made 
an account of it on the basis of the envelope theory [32] . 
Eq. (2.9) leads to the equations for A(t) and </>(t), 



A = e- 1 - 

2 V 4 



(2.10) 



which equations may be called an amplitude and phase 
equation, respectively. The original equation is reduced 
to these simpler equations for the amplitude and phase 
which parametrize the solution of the original equation. 
These reduced equations are readily solved, with which a 
resummation of the perturbation series is performed; the 
resumed solution is found to successfully admit a limit 
cycle with a radius of 2. 

The resultant envelope function as a global solution is 
given by 

x E (t) = x(t;t) = W(t) 

= A{t) cos(j)(t) - e4J^sin30(t), (2.11) 

with A(t) and <f>{t) being the solution of Eq. (2.10). Thus, 
we have succeeded in not only obtaining the asymptotic 
solution as a whole but also extracting the slow variables 
A{t) and 4>(t) explicitly and their governing equations. 

However, there is a problem left: Does x E (t) = x(t;t) 
indeed satisfy the original differential equation (2.1)? We 
give here a proof for that [32] . 

First let us rewrite Eq. (2.1) into a coupled equation 
of first order. 



dq(t) 
dt 



where 



Q 
F 



*((7i = x, q 2 = x), 



(2.12) 



(2.13) 
(2.14) 



We have an approximate local solution to Eq. (2.12) q(t; t ) 
around t — t up to 0(e n ), corresponding to x(t;t ); 



^ = F( q (t;to);e) + 0(e n ). 

Now, the RG/envelope equation implies that 
dq(t;t ) 



dt 



t=t 



= 0. 



(2.15) 



(2.16) 



The envelope function q E (t) corresponding to x E (t) is de- 
fined by 



<7eM = q(t-,t). 



(2.17) 



It is now easy to show that q E (t) satisfies Eq. (2.12) up 
to the same order as q(t; to) does: In fact, Vt = to 



dq E (t) 



dt 



dq(t;t ) 



t=t a 



dt 



dq(t;t ) 



+ 



dq(t;t ) 



t=t 



dt 



t=t 



dt 

t=t a 

F(q(t;t);e) + 0(e n ) 
F(q E (t);e)+0(e n ). 



(2.18) 



This completes the proof. Here Eq. (2.16) and Eq. (2.15) 
have been used together with the definition of q E (t), Eq. 
(2.17). It should be stressed that Eq. (2.18) is valid uni- 
formly Vt i.e., in the global domain of t, in contrast to Eq. 
(2.15) which is in a local domain around t — to- 

We can summarize what we have done as follows: when 
there exist zero modes of the unperturbed operator, the 
higher-order corrections may give rise to secular terms, 
which are renormalized into the integral constants in the 
zeroth-order solution using the RG/envelope equation, and 
thereby, the would-be integral constants are lifted to dy- 
namical variables. 

It will be found that the would-be integral constants 
corresponding to A and <fi exactly correspond to the hy- 
drodynamic variables characterizing the local equilibrium 
distribution function, such as the temperature T, chemical 
potential fi, and flow velocity (w M w M = 1): The equa- 
tions (2.10) of the amplitude and phase which parametrize 
the solution also exactly correspond to the hydrodynamic 
equation governing the hydrodynamic variables which 
parametrize the distribution function as the solution of 
the Boltzmann equation. 



3 Relativistic Boltzmann equation 

The relativistic Boltzmann equation (RBE) reads [20,21] 

P^f p (x) = C[f] p (x), (3.1) 

where f p (x) denotes the one-particle distribution function 
with jj m being the four-momentum of the on-shell particle, 
i.e.,p^Pft = p 2 = m 2 andp > 0. C[f] p (x) inr.h.s. denotes 
the collision integral 



,Pl\P2,P3) 



c[f] P (x)= ^EiEiEi^ 

PI 1 P2 2 P3 3 

x (fp2 (x)f P3 (&) - / P (aO/pi (»)) , (3.2) 



where ui(p , p\\pi , ^3) denotes the transition probability 
due to the microscopic two-particle interaction with the 
symmetry property 

w(p,Pl|i>2,P3) = w (P2,P3b,Pl) 

= w(pi,p|p 3 ,p 2 ) = w(p 3 ,p 2 |pi,p),(3.3) 
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and the energy-momentum conservation 

Pi\P2 , Ps) oc <5 4 (p + pi -p 2 -ps). 



(3.4) 



To make explicit the correspondence to the general formu- 
lation of the reduction theory given in [35,33], we treat the 
momentum as a discrete variable; apart from such a formal 
reasoning, the summation with respect to the momentum 
may be interpreted as the integration in practical use as 
follows, 



E - fa 



(3.5) 



with q being the spatial components of the four momen- 
tum q 11 . 

For an arbitrary vector (p p (x), the collision operator 
satisfies the following identity thanks to the above- 
mentioned symmetry property, 

^2^o<Pp(x)C[f\ p (x) 
p 

= E O 1 ^(P)Pl|P2,P3) 

2! 4 p^pipipi 

P,Pl~P3 

x (<p p (a;) + (/? Pl (a;) - tp P2 (x) - <£ P3 (x)) 

x (/p2 (a:) /p 3 (») ~ f P ( x ) fpi (a:)) • (3-6) 

Substituting (1, p M ) into <p p (a;) in Eq. (3.6), we find that 
(1, p^) are collision invariants satisfying 



due to the particle-number and energy-momentum conser- 
vation in the collision process, respectively. We note that 
the function ipo p (x) = a(x) + p^ b fl (x) is also a collision 
invariant where a(x) and b fi (x) are arbitrary functions of 
x. This form is, in fact, the most general form of a collision 
invariant [20]; see [21] for a proof. 

Owing to the particle-number and energy- 
momentum conservation in the collision process leading to 
Eq. (3.7), we have the balance equations for the particle 
current N )1 (x) and the energy-momentum tensor T fiu (x) 7 



d^N^x) EE 0, 



d v T^{x)=d v 



E i^/pO*) 



0, (3.8) 



(3.9) 



respectively. It should be noted that any dynamical prop- 
erties are not contained in these equations unless the evo- 
lution of f p (x) has been obtained as a solution to Eq. 
(3.1). 

In the Boltzmann theory, the entropy current may be 
defined [20] by 



S " (*) = ~ E [in ((2 f P (x)) 1 



where the factor (27r) 3 is necessary owing to our conven- 
tion (3.5) [20]. The entropy current S^(x) satisfies 

d^(x) = -J2^C[f} p (x) ln((2^) 3 / p (x)),(3.11) 



due to Eq. (3.1). One sees that S^(x) is conserved only if 
ln((2 7r) 3 f P (x)) is a collision invariant, i.e., ln((2 7r) 3 f p (x)) = 
l -Po P {x) — a(x) + p^ b^(x). One thus finds [20,21] that 
entropy-conserving distribution function may be parametrized 



as 



f P (x) 



(2tt)* 



exp 



T(x) 



f e p q (x), 



(3.12) 



with u^(x)u fl (x) — 1. The function (3.12) is identified 
with the local equilibrium distribution function called the 
Jiittner function [43], where T(x), fi(x), and u^(x) in Eq. 
(3.12) are the local temperature, chemical potential, and 
flow velocity, respectively; see [21] for a proof. These five 
variables are called hydrodynamic variables. Due to the 
energy-momentum conservation in the collision process, 
we see that the collision integral identically vanishes for 
the local equilibrium distribution f p q (x) as 



C[f e \(x) = 0. 



(3.13) 



Some remarks are in order here. In the proof [21], the 
Gibbs-Duhem relation as given by 



d(s/n) = ^ (d(e/n) 



(3.14) 



(3.10) 



is taken for granted, where s, e, n, and p denote the en- 
tropy, internal energy per volume, particle density, and 
pressure in the equilibrium state, respectively. However, 
Van and Biro [17] have recently argued that the conven- 
tional Gibbs-Duhem relation (3.14) may be modified so 
as to contain the contribution from the thermal flow in 
the local equilibrium state of a relativistic system, and 
given a different interpretation for T(x), p,{x), and u^{x) 
in (3.12); this modified definition of the local equilibrium 
state, they claim, leads to the relativistic hydrodynamic 
equation in the particle frame with the stable equilibrium 
state. Although this is certainly an interesting possibility, 
we will not follow this novel interpretation in this review: 
We shall make some comments on some related problem 
below. 



4 Reduction to hydrodynamic equation 

Let us try to solve the RBE (3.1) in the hydrodynamic 
regime, and thereby derive the hydrodynamic equations 
governing the hydrodynamic variables. 
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4.1 Relativistic Boltzmann equation in local rest 
frame of flow velocity 

To make it explicit to solve the RBE in the hydrodynamic 
regime, we first convert the RBE (3.1) into the following 
form with the use of the flow velocity (u^u^ = 1) [20] : 

|l/p(r, a) = ^C[f] p (r, a) 



1 



p ■ u 



p-V/ p (r, a), (4.1) 



where the new coordinate system (t, er M ) is defined as fol- 
lows, 



d 



(4.2) 



= ( g ^ - u ^ u ^)d u =A^d v = y^. (4.3) 

We note that D and arc temporal and spatial differ- 
ential operators familiar in the literature. In Eq. (4.1), the 
small parameter e is introduced as a measure of the non- 
uniformity of the fluid, which may be identified with the 
Knudsen number; e will be set back to unity in the final 
stage of the analysis. In the present analysis based on the 
RG method, the pcrturbative expansion of the distribu- 
tion function with respect to e is first performed with the 
zcroth-order being the local equilibrium one; the dissipa- 
tive effect is taken into account as a deformation of the 
distribution function made by the spatial inhomogeneity 
as the perturbation. Thus the above rewrite of the equa- 
tion with e reflects a physical assumption that only the 
spatial inhomogeneity is the origin of the dissipation. It 
is noteworthy that our RG method applied to the nonrel- 
ativistic Boltzmann equation with the corresponding as- 
sumption successfully leads to the Navier-Stokes equation 
[45]; the present approach [24,26] is simply a relativistic 
generalization of the nonrelativistic case. 

Here we make a comment on the possibility of a rewrite 
of the RBE (3.1) with use of a different time-like four 
vector in place of the flow velocity w M . In other words, we 
examine whether Eq. (4.1) with u M being identified with 
the flow velocity is a unique rewrite of the RBE (3.1) in 
a covariant manner. We argue that it is the case on the 
basis of a physical ground. 

We first introduce a generic time-like four vector a M 
with a 2 > 0, and call it the macroscopic- frame vector, 
following [24, 37] . Without a loss of generality, the generic 
vector of Lorentz covariance takes the form, 

= A 1 u" + A 2 d»T + A 3 d^fj, + A A u v d v u», (4.4) 

since and 9 M are the only available Lorentz vectors at 
hand. Here, A\, A 2 , A 3 , and A 4 arc arbitrary functions 
of the temperature T and the chemical potential /i; Ai = 
Ai(T, fj) for i = 1, 2, 3, 4. Owing to the identity 

d» = u^u v d v + {g^ - u^u v )d u = u^D + V, (4.5) 



where D and V M have been defined in Eq.'s (4.2) and 
(4.3), Eq.(4.4) is rewritten as 

= {A l + A 2 DT + A 3 Dfj) 

+ A 2 V"T + A 3 VV + A 4 Du» 
= C t (T, n) u» + A 2 V^T + A 3 VV + A 4 Du^, 

(4.6) 

with C t (T, (i) = A 1 + A 2 DT + A 3 D\l. The relative mag- 
nitudes of Ct and ^2,3,4 are only constrained by the in- 
equality a 2 > in the present stage. However, it should be 
emphasized that the space-like terms with the coefficients 
^2,3,4 are all derivative terms, which are supposed to be 
small in the hydrodynamic regime even in the dissipative 
regime if the dynamics is governed the hydrodynamics at 
all. 

By replacing by a M given by (4.6), we have the 
generic coordinate system (f , a^) as defined by 



d c 

OT 1 



(4.7) 
(4.8) 



With this coordinate system, the RBE is rewritten as 



d 



1 



- e F^fl§^< f '*>' (4 ' 9) 

where the e is again multiplied to <9/<9cr M as was done 
in Eq. (4.1) where it is supposed that only the spatial 
inhomogeneity is the origin of the dissipation. Then the 
space-like terms with the coefficients A 2 and A 3 in a M are 
of higher order with respect to e and should be ignored 
in this set up. Furthermore, since we start with a station- 
ary solution with vanishing time-dependence in the RG 
approach, the term with A 4 should be also ignored. Thus 
we have 



a» = C t (T, n)u» = bu^, 



and accordingly 



d 



= V M . 



(4.10) 



(4.11) 



If we naturally require that b should be independent of the 
momentum p^ , it is easy to show [44] that the "normaliza- 
tion" factor b can be made unity without loss of generality, 
in conformity of the natural choice [20,21] = u^. 

It is remarkable that this natural choice uniquely leads 
to the hydrodynamic equation in the energy (Landau- 
Lifshitz) frame, as will be shown and discussed later [24, 
37]. Conversely, a choice of b different from unity with 
a momentum dependence could lead to various hydrody- 
namic equations other than that of Landau and Lifshitz, 
including the one in the particle frame for viscous flu- 
ids as was shown by the present authors [24,37]. How- 
ever, it is worth emphasizing that the particle frame can 
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be only realized when b has a peculiar momentum de- 
pendence such as b = m/(p ■ u) (a M = {m/(p ■ u))u^) 
[24,37]. In retrospect, however, the possible momentum 
dependence of b can not be legitimate for a M to play a 
macroscopic-frame vector, because it means that the co- 
variant and macroscopic space-time in the particle frame 
is defined for a respective particle state with a definite 
energy-momentum, which is certainly unnatural and lead 
to a trouble in a physical interpretation [44]. Thus, we 
naturally require that b is independent of the momentum 
p^ and hence a M — u^. 



4.2 Hydrodynamics from relativistic Boltzmann 
equation by renormalization-group method 

Applying the perturbation theory to Eq. (4.1), we de- 
rive the relativistic dissipative hydrodynamic equation as 
the infrared asymptotic dynamics of the RBE by the RG 
method [32,33,24,26]. 

In this approach, we first try to obtain the perturbative 
solution f p to Eq. (4.1) around the arbitrary initial time 
t = tq with the initial value / p (tq, a); 



Ip(t = t ,<t;t ) = f p (T ,a). 



(4.12) 



Note that the solution depends on the initial time To at 
which f p (r = t , cr; t ) is supposed to be on an exact so- 
lution. We expand the initial value as well as the solution 
with respect to e as follows: 

/ p (t,o-;t ) = ./^ 0) (t,ct;to) +£/^ 1) (t,o-;t ) 

+£ 2 /^(t, ( t;to) + ---, (4.13) 

The zeroth-order equation reads 
1 



(4.14) 



C[/(°)] p (T,a;T ). (4.15) 



p ■ u 



Since we are looking for the slow motion to be realized 
asymptotically when t — > oo, we take the stationary solu- 
tion satisfying 



^/(°)( T , CT ; T0 ) =0 , 



(4.16) 



implying that C[f^] p (r, a; To) = Vcr, which is solved by 
a local equilibrium distribution function, i.e., the Jiittner 
distribution function, 



/(°)(T,a;T ) = 



exp 



(2tt) 3 



Mo-; t ) - p 11 ^(ctjto) 



T{<j;t ) 



(4.17) 



with u^(<r; tq) u M (cr; To) = 1. Here the would-be integra- 
tion constants T(a ; To), ^(cr ; t ), and u ll (a ; t ) are inde- 
pendent of t but may depend on to as well as cr. 



Now that the zero-th order solution is given, the first- 
order equation reads 

^\r) = Y,ApJ^\r)+Fp, (4.18) 



with 



Fp = - 



p ■ u 



cq 



(4.19) 



where A pq denotes a matrix element of the linearized col- 



lision operator A; i.e., 



(A)pq — A pq 



1 



d 



p-u df q 



C[f] f 



(4.20) 



Let us examine the spectral properties of A; for which, 
it is found convenient to convert A to another linear op- 
erator, 



L^iD^An, 



(4.21) 



with the diagonal matrix (f cq ) pq = ,f p q S pq . Next we define 
an inner product between arbitrary nonzero vectors ip and 
ip by 



, 4)(p- u )/p>pV> p , (4-22) 



which satisfies the positivc-dcfinitencss of the norm as 

^^=Ei(p^)/p c, w 2 >° ( 4 - 23 ) 



for ifp 0, since both p^ and w M are time-like vectors with 
p° > 0. 

Then it can be shown [24,26] that the linearized colli- 
sion operator L has remarkable properties that it is semi- 
negative definite and has five zero modes given by 



for a = fi, 



1 x in for a = 4. 



(4.24) 



The functional subspace spanned by the five zero modes 
is called the Po space and the projection operator to it is 
denoted by Po! 

[P^]p = <Pz P %^{4^), (4-25) 

where %^ is the inverse matrix of the the Po-space metric 
matrix rfc 13 defined by 



(4.26) 



We also call the complement to Po the Q space and in- 
troduce Qo = 1 — Po- In the following, we also use the 
modified projection operators defined by 

Po - f cq Po(.rr\ Qo = / eq Qo(/ eq )-\ (4.27) 



Kyosuke Tsumura, Teiji Kunihiro: Title Suppressed Due to Excessive Length 



7 



which means, for example, 

[Po V] p = fFv&V^pirf > (/ eq )-V). (4-28) 

Then the perturbative solution up to the second order 
reads 



/p(t,<t;t ) 



/W(T,a;ro)+e/«(r, a;r ) 
+ £ 2 /f(r, ( 7;To) + 0(£ 3 ), (4.29) 



where / (1) (r, <r; r ) = (t-t )Po^-4 _1 Qo-F and a lengthy 
formula for f^ 2 \r,a;To), which we do not write down for 
the sake of space; see [26] for the details. 

We remark that this solution contains secular terms, 
which apparently invalidates the perturbative expansion 
for t away from the initial time To- We can, however, uti- 
lize the secular terms to obtain an asymptotic solution 
valid in a global domain [32,33]. Indeed we have a fam- 
ily of curves / p (r, cr; r ) parameterized with r : They are 
all on the exact solution f p (a ; r) at r = To up to 0(e 3 ), 
although only valid for r near To locally. Then, the enve- 
lope curve of the family of curves, which is in contact with 
each local solution at r = t , will give a global solution in 
our asymptotic situation, which is shown to be the case 
[32,33]. According to the classical theory of envelopes, the 
envelope that is in contact with any curve in the family 
at t = r is obtained [32] by 



d ~. 
dr 



fp{T,a;T ) 



= 0. 



(4.30) 



The derivative w.r.t. To hits the hydrodynamic variables, 
and hence we have the evolution equation of them that 
is identified with the hydrodynamic equation [45,24]. We 
also note that the invariant manifold which corresponds to 
the hydrodynamics in the functional space of the distribu- 
tion function is explicitly obtained as an envelope function 
[24,26]: /e p (t,ct) = / p (t, a ; To = r), the explicit form of 
which is referred to [24,26]. We note that this solution is 
valid in a global domain of time in the asymptotic region 
[26]. 

Putting back e to 1, Eq. (4.30) is reduced to the fol- 
lowing form in this approximation, 



(L eq + <5/«) = 0.(4.31) 



where Sfp denotes the first-order correction to the dis- 
tribution function 



5fW = -[A^QoFlp. 



(4.32) 



If one uses the identity (p ■ u) d/dr + p ■ V = p^ d^, Eq. 
(4.31) is found to have the following form 



with = + 6T& and N? st = M )" + S N? st . 

Here, 

T (o)^ = 4rPV/ p cq = e// -pA» v , (4.34) 
p 

AT(°^ = ^l^/7 = n^, (4.35) 

p 

with e, p, and n being the internal energy, pressure, and 
particle-number density for the relativistic ideal gas, re- 
spectively, while the dissipative parts are given as a devi- 
ation of the local equilibrium distribution function 

^^E^wi 1 '' ( 4 - 36 ) 
p 

57V lst =^l^/«. (4.37) 
p P 

As is well known, the local equilibrium distribution func- 
tion as given by (4.17) only gives the relativistic Euler 
equation without dissipation. 

4.3 Possible uniqueness of Landau-Lifshitz frame 

In this subsection, we present the explicit form of the dis- 
sipative parts ST^ t and SN^ st and discuss their properties. 
An evaluation of Eq. (4.36) together with (4.32) gives [24, 
26], 



5Tj% = C A» v V • u + 2 n A^ pa V 



h? T 



(4.38) 
(4.39) 



respectively, with A» v p° = 1/2 • (A^A V(y + A pa ' A v ? - 
2/3 • A pv A 9 "). Here, h denotes the reduced enthalpy per 
particle. The bulk and shear viscosities and the thermal 
conductivity are denoted by £, n and A, respectively. It 
is clear that these formulas completely agree with those 
proposed by Landau and Lifshitz [16]. Indeed, the respec- 
tive dissipative parts 8T^ t and 5N^ st in Eq.'s (4.38) and 
(4.39) meet Landau and Lifshitz's ansatz 



Se = <5Tf s j u v = 0, 
5n = u^dN^ t = 0, 
Q» = A^ v 6T"£u p = 0. 



(4.40) 
(4.41) 
(4.42) 



Thus we find that the frame on which the flow velocity is 
defined inevitably becomes the Landau-Lifshitz (energy) 
frame, if the hydrodynamics is to be consistent with the 
underlying relativistic Boltzmann equation . 

Let us see the above fact in the level of the distribution 
function. We first note that Eq. (4.36) and (4.37) can be 



dpTZ = 0, « t = 0. 



(4.33) 



1 The uniqueness of the energy frame for the relativistic hy- 
drodynamics is recently argued also in a different context [46] . 
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rewritten as 



with 



= E i>Mf?*p> 

V 

^ = Ei^W< 



(4.43) 
(4.44) 

(4.45) 



which belongs to the Q space and thus orthogonal to the 
zero modes, 



(<p% , </>) = for a = 0, 1, 2, 3, 4. 



(4.46) 



Here, the inner product is defined by Eq. (4.22). Then, 
Eq. (4.46) with a = \i is reduced to 







= u v 5TZ- 



W )/ p e v^ = ^Ei^W 



(4.47) 



Similarly, Eq. (4.46) with a = 4 is reduced to w M SN^ st = 
0. Thus, one can readily see that these equations coincide 
with Landau and Lifshitz's ansatz: We remark that Eq. 
(4.47) implies the following two equations, Se = u v 5T^ t 
= and Q p = A pp u v 5T^ t = 0, which are nothing but 
the matching conditions [20] imposed to select the energy 
frame in all the other existing approaches based on the 
Boltzmann equation. In other words, we have given the 
foundation to the matching conditions [20] for the energy 
frame. 

We can present an intuitive picture of why the energy 
frame is uniquely selected in this method. As shown in 
[37], the physical quantity transported by each particle 
governed by the RBE (4.1) can be identified with (p ■ u). 
To clarify the physical meaning of (p ■ u), we take the 
non-relativistic limit of this quantity: 



(p ■ u) ~ m + 



rn 



P_ 

P° 



— u 



(4.48) 



where — (u°, u) andp M = (p°, p). This equation shows 
that {p ■ u) can be interpreted as the kinetic energy of the 
fluid component measured in the rest frame of u^. Thus, 
it is natural that the resultant equation becomes the one 
in the energy frame adopted by Landau and Lifshitz. 

A remark is in order here. The uniqueness of the en- 
ergy frame comes from the two natural conditions used in 
the derivation, i.e., the identification of the time- like vec- 
tor in the Jiittner distribution function (4.17) with the 
flow velocity and the physical assumption that the dissi- 
pativc effect comes from only the spatial inhomogeneity. If 
one of these conditions were to be challenged, as claimed 
in [17], for instance, the uniqueness of the energy frame 
could be violated. It is clear that further studies are needed 
for establishing the uniqueness of the energy frame in the 
relativistic hydrodynamics for viscous fluids. 



4.4 Transport coefficients 

Since our theory starts from a microscopic theory as sta- 
tistical mechanics, we have the microscopic expressions for 
the transport coefficients appearing in the hydrodynamic 
tensor(4.38) and current (4.39), as follows: 



C=-^(fl,L- 1 II), 

A = ^(j",r 1 J (i ), 



(4.49) 
(4.50) 
(4.51) 



Here, we have introduced the following microscopic ther- 
mal forces (it p , Jg, 7f£") ee (7T p , J£, 7r£")/(p • «), with 

n p = (4/3 - 7) (p • u) 2 + (( 7 - l)Th - 7 T) (p ■ u) 

-1/3-m 2 , (4.52) 
Jg = -((p.u)-Th)A^p„, (4.53) 
tt^ = A^p pPa . (4.54) 

Here, 7 denotes the ratio of the constant pressure and 
volume heat capacities. We note that the microscopic ex- 
pressions for the transport coefficients (4.49)-(4.51) are 
in agreement with those given by the Chapman-Enskog 
method [20]. 

It is noteworthy that the transport coefficients can be 
rewritten in the Green-Kubo formula [28]. With the use 
of the "time-dependent thermal force" defined by 



n p (s) = J2[e°% q n q 



(4.55) 



and so on, the relaxation functions are given by the time- 
correlators 



fl f ( s ) = 1(77(0), J7(*)>, 



(4.56) 



and so on for R\(s) and R n (s) with obvious modifications 
to Rq(s). Then the transport coefficients given in Eq.'s 
(4.49)-(4.51) are rewritten as follows [24,26], 

/>oc />oo />oo 

C = / dsR^(s), A= / dsR\(s), rj = dsR n (s). 
Jo Jo Jo 

(4.57) 



5 Generic stability of relativistic 
hydrodynamic equation in energy frame 

In this section, we shall provide a proof [26] that generic 
constant solutions of the relativistic dissipative hydrody- 
namic equation in the energy frame is stable against a 
small perturbation [47], on account of the positive defi- 
niteness of the inner product as shown in Eq. (4.23). 
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For this purpose, we first note that F p is reduced to 



Fn 



-f 

J 1 



eq 



p ■ U 



(5.1) 



with <£q p being the zero modes defined in (4.24) and 

{—u u /T for a = v, 
m' 1 fi/Tioi a = 4. 
Then Eq. (4.31) is rewritten in the following form, 



(5.2) 



./p oq (l + 



[L-^f^Xp) 



where 



<f%=[Q <fff p , 



(5.3) 



(5.4) 



with iflp =p^tp^ p /(p-u). 

Now, a generic constant solution means that it de- 
scribes a system having a finite homogeneous flow with a 
constant temperature and a constant chemical potential, 
as follows: 

T(ct;t) = T , ^(ct;t) = fj, , u M (cr;r) = u p, (5.5) 

where To, fiQ, and uq^ are constant. We note that these 
states include the thermal equilibrium state as a special 
case. 

We shall show the linear stability of the constant solu- 
tion of the relativistic dissipative hydrodynamic equation 
in the energy frame. We represent T, fj,, and around 
the constant solution as follows: 



T(a; r) = T + ST(a; r), 
fj,(a; t) = fiQ + 6fj,(a; t), 
u^(<t ; t) = it 0M + Su^a ; t), 



(5.6) 
(5.7) 
(5.8) 



where the deviations ST, 5fi, and <5w M are assumed to so 
small that terms in the second or higher orders of them 
can be neglected. Instead of these six variables which not 
independent of each other because 5u^ u 1 ^ = 0, we use the 
following five independent variables, 



SX n 



(-S(^)=- S -^+ST^ fora = , 
\T J T T n 2 P ' 



rn 



m- 1 & - ST* 
K T a T 2 



) for a = 4. 



(5.9) 

Substituting Eq. (5.9) into Eq. (5.3) and with some ma- 
nipulation, we obtain the linearized equation governing 
SX a as 

({^,4) + {y«,L-^f)V v )^5Xe 



+ 



(5.10) 



Here, we have used the following simple relation 

6(fZ q ) = f c p \o p SX a , (5.11) 

We note that all of the coefficients in Eq. (5.10) take a 
value of the constant solution (T, fi, u M ) = (To, no, uq^). 
Owing to the orthogonality between the Po and Qo spaces, 
Eq. (5.10) is reduced to 

vf ^SX + B^SX^O. (5.12) 

Here ?7q^ is the metric tensor defined in (4.26) and B a P is 
defined by 

^^r^v^r^v,, (5.i3) 



with 



rfr* = «\£-vr>. (5.i4) 

Both 770 and B are symmetric tensors. 

With the ansatz '0^(0-57-) = 5X a {k; A)c tk - a - AT , Eq. 
(5.12) leads to the following algebraic equation, 



{Ar,f-B a ")5Xp=0, 



(5.15) 



with B ali = i{<Pi a , <Pq)/c m — rft av ^k tl k v . Thus we have the 
eigenvalue equation as follows, 



det(yl?7o - B) = 0, 



(5.16) 



which would give the dispersion relation A = A(k). The 
stability of the generic constant solution (5.5) against a 
small perturbation is assured when the real part of A(k) 
is nonnegative for any fc M , which is shown to be the case 
as follows. 

Now, recall that the metric matrix 770 is a real sym- 
metric and positive-definite matrix, which implies that it 
has a Cholesky decomposition, 

! -*UU, 



Vo 



(5.17) 



where U denotes a real matrix and U a transposed matrix 
of U. Then Eq. (5.16) is converted to 



dct(AI -U B'U) =0, 



(5.18) 



where I denotes the unit matrix. Eq. (5.18) tells us that 
A(k) is an eigen value of U B l U. 

There is a following theorem: The real part of the eigen 
value of a complex matrix C is nonnegative when the Her- 
mite matrix Rc(C) = (C + C')/2 is semi-positive defi- 
nite. Applying this theorem to the present case, we find 
that the real part of A(k) becomes nonnegative for any /c M 
when Re(U B t U) is a semi-positive definite matrix, which 
is shown to be the follows; 

= w a [URc(BYU] ap w 

= [wU] a [Re(B)] al3 [wU]p 

= -[wU] a r,r^ K k v [w U] p 

= -( fc M [w U] a (pf a , L^ 1 k v [w U]f) (fii 13 ) 

= -(V, £ _1 V>) > for w a + 0, (5.19) 
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with ip p = k^[wU] a (fip . This completes the proof that 
the generic constant solution in Eq. (5.5) is stable against 
a small perturbation. 



6 Second-order equations and moment 
method 

In the first-order hydrodynamic equations, the zero modes 
of the linearized collision operator form the invariant man- 
ifold on which hydrodynamics is defined; the would-be 
constant zero modes acquire the time-dependence on the 
manifold by the RG equation. Our formalism can be ex- 
tended so as to incorporate excited modes as additional 
components of the invariant/attractive manifold [36], and 
hence we can derive an extended thermodynamics or Israel- 
Stewart type equation with novel microscopic expressions 
of the relaxation times and lengths [36] . Furthermore, our 
theory suggests a proper ansatz for the distribution func- 
tion to be used in the moment method [36] . For the short- 
age of space, we here give a sketch of some of our results for 
the extended thermodynamics, leaving the detailed deriva- 
tion in a separate paper [36]. We emphasize that our the- 
ory gives an explicit construction of the invariant manifold 
corresponding to thirteen moments, which has been long 
sought for [39,40]. 



6.1 A brief review of Grad's thirteen-moment method 
and Grad-Muller equation: non-relativistic case 

In Grad's thirteen-moment method [42,20], the one-particle 
distribution function fv(t, x) is represented as 



Thus, the Grad-Muller equation is obtained as a closed 
system of the equations governing T, n, u l , 7r* J , and J 1 
in terms of the transport coefficients and the relaxation 
times, which are given in terms of an inner product defined 
by (ip, X) cq = J2v f^^vXv- For example, the shear 
viscosity and the relaxation time of the stress tensor n t: > 
are expressed as 



f v (t,x) = f^(t,x) (! + #„(*, as)), 



(6.20) 



where denotes the Maxwell distribution function and 
$v the deviation from /£ q given by 

$v(t,x) = n%(t,x)Tr ij (t,x) + J l v {t,x)J l {t,x) 

= ^(t,x), (6.21) 

with 

7rj/(t, x)=m (5v% x)Sv j (t, x) -^S ij \Sv(t, x)\ 2 ^j, 

(6.22) 

4(t, x) ee \6v{t, *)\ 2 -\ T{t, *)) 8v% x). (6.23) 

Here, Sv(t, x) = v — u(t, x) is the peculiar velocity. 

Then the evolution equation of the thirteen coefficients 
are determined by the equations all of which are derived 
from the Boltzmann equation (d/dt + v ■ V)/y(t, x) = 
C[f]v(t, x) with use of the linearized collision operator 
given by 



j-eq 



(6.24) 



1 (TC ij ,TT ij ) 



7T 

cq \ ' " /eq 



10T (tt 



cq 



[7T kl , L% M ) 



(6.25) 
(6.26) 



It is well known that the formula (6.25) is different from 
that given in the Chapman-Enskog expansion method, 
and there are many attempts both in non-relativistic [40] 
and relativistic cases [23,48] to modify and/or extend the 
Grad moment method so that the transport coefficients 
thus obtained become consistent with those obtained by 
Chapman-Enskog method. 



6.2 Relativistic mesoscopic dynamics from the RG 
method 

In this subsection, we make a brief report on our attempt 
[36] to extend the RG method so as to obtain the so called 
mesoscopic dynamics [39] in the relativistic case. 

First we show the results in such a way that a com- 
parison with the Grad moment method is apparent. If we 
express the distribution function by f p (x) — fp q (x) (l + 
<&p(x)), our RG method gives the following expression of 



- VL 

y PI 
<i 



(6.27) 



f=F 



where 7T p , Jp 1 , and tv£ v are the microscopic thermal forces. 
This new form is different from any proposals in the liter- 
ature [23,48]. 

The resultant energy-momentum tensor and particle 
current are found [36] to have the following forms, respec- 
tively, 

T 2nd = e U ^ u ' y -(P + n ) Ailv + . (6-28) 
K nd = nu» + J». (6.29) 

The relaxation equations derived in our RG method read 

77 = -CV -u-rnDII 

+ other terms involving relaxation lengths, (6.30) 

J» = \±-V^-rjA» a DJ a 
h A 1 

+ other terms involving relaxation lengths, (6.31) 
Tr"" = 2 77 A^ pa V p u a - TV A^ ab Dn ab 
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+ V • U + K £) A hpa AaMe yd u e 



1 per 



+ other terms involving relaxation lengths, (6.32) 

where uj^ = \ ( W - W) is the vorticity. 

Our RG method [36] gives microscopic expressions for 
the relaxation times 777, tj, and iv as follows 



(77, L 

Tn = - 



•■n) _ j °°d s s r c (s) 



(n,L-^n) J^ds r c ( s ) 



(6.33) 



and so on for tj and r T with obvious modifications. We 
note that our novel formulae for the relaxation times are 
all nicely represented in terms of the relaxation functions 
R^(s), R\(s), and i?,,(s) so that they have a natural phys- 
ical meaning of the relaxation time as the correlated time 
of the respective relaxation function in contrast to other 
approaches [23,42,48]. We also mention that bulk and 
shear viscosities and the heat conductivity derived in our 
method do coincide with those in the Chapman-Enskog 
method as shown before. 



7 Summary and concluding remarks 

We have reported our attempts to derive first-order and 
second-order relativistic hydrodynamic equations from rel- 
ativistic Boltzmann equation which has a manifest Lorentz 
invariance and does not show any pathological behavior 
such as the instability and acausality seen in existing hy- 
drodynamic equations. We have given an argument, on a 
physical ground on the nature of the origin of the dissi- 
pation and the form of the local equilibrium distribution 
function, that the energy frame is uniquely chosen as the 
one in which the relativistic hydrodynamic equation for a 
viscous fluid is defined. We have given the novel extended 
thermodynamics both in non-relativistic and relativistic 
cases through the explicit construction of attractive man- 
ifold containing the relaxation process from Boltzmann 
equation. 

It is worth emphasizing that that all the equations 
derived in this work are consistent with the underlying ki- 
netic equation, i. e. , relativistic Boltzmann equation. This 
is one of the advantage in our theory because our the- 
ory explicitly gives the solution (distribution function) 
of the Boltzmann equation, which is expressed with the 
hydrodynamic variable and relaxation times and lengths, 
and thereby makes a systematic description of the time- 
evolution of the system from hydrodynamic to kinetic 
regime. Such an overall analysis should be desirable for 
that of the matter created at RHIC, LHC and other sys- 
tems where the proper dynamics would change from the 
hydrodynamic to the kinetic ones or vice versa [4-9] . Fur- 
thermore, it would be interesting to evaluate the relax- 
ation times as well as the transport coefficients of the cre- 
ated matter with the use of the microscopic representa- 
tions obtained in this work. 

Finally, we note that the renormalization-group method 
[31-33,45] itself has a universal nature and can be applied 



to derive a slow dynamics from kinetic equations other 
than the simple Boltzmann equation, say, Kadanoff-Baym 
equation [49]. 
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